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MAGNETIC FIELD 
By Vernon J. Rossow 


SUMMARY 


The effectiveness of a magnetic field in stabilizing the laminar 
flow of an incompressible, electrically conducting fluid is studied. 

The neutral stability curves pertaining to a two-dimensional sinusoidal 
disturbance are presented for flow over a semi- infinite flat plate in 
the presence of either a coplanar or transverse magnetic field and for 
channel flow in the presence of a coplanar magnetic field. As is to be 
expected, the magnetic field stabilizes the flow unless the velocity 
profile is distorted by the magnetic field to an inherently unstable 
shape. This occurs when a transverse magnetic field is fixed relative 
to a semi- infinite flat plate. 


INTRODUCTION 


Mere mention of the possibility of controlling the motion of 
electrically conducting fluids with a magnetic field stimulates one’s 
imagination to conceive flow fields which may furnish certain ideal char- 
acteristics. All too often the configurations are too complicated to be 
amenable to analysis and one must be content with a greatly simplified 
version of the original idea. A survey of the literature shows that a 
number of basic solutions are being accumulated. A large portion of the 
effort is directed at the theoretical evaluation of the effectiveness of 
a magnetic field in stabilizing a given laminar flow so that transition 
to turbulent flow is inhibited. Some of the earliest work on problems 
of this type was carried out by S. Chandrasekhar. He found that a mag- 
netic field would inhibit the onset of convection in a fluid heated from 
below (ref. 1), and would impede the transition to turbulence of fluid 
between rotating cylinders of nearly the same diameter (ref. 2). In a 
later paper, reference 3> it is found that a layer of fluid heated from 
below and subject to rotation is, under certain conditions, destabilized 
by application of a small magnetic field. The motion is stabilized by 
increasing the magnetic field strength beyond a certain amount. 

The effect of a magnetic field on the stability of the flow of an 
incompressible electrically conducting fluid in a two-dimensional channel 
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has been studied for a coplanar magnetic field by Stuart (ref. 4) and for 
a transverse magnetic field by Lock (ref. 5). The transverse magnetic 
field is found to be the more effective in stabilizing the flow field. 

The high degree of stabilization brought about is attributed (to the 
order of accuracy of the analysis) entirely to the change in the velocity 
profile caused by the interaction of the fluid and magnetic field. When 
the magnetic-field lines are parallel to the stream direction, the favor- 
able effect on the stability of a disturbance is brought about by the 
electromotive resistance encountered when a fluid element leaves its 
normal path of motion in an effort to form a turbulent eddy, thereby 
crossing magnetic lines of force. 

The effect of a coplanar magnetic field on the stability of a laminar 
mixing region was studied by Curie (ref. 6). The Reynolds numbers at 
which a small disturbance becomes unstable are generally quite small for 
this type flow field (generally less than 100) but increase rapidly with 
increasing magnetic parameter. Complete stabilization is predicted for 
a magnetic parameter over 0.301. 

An experimental example of flow instability caused by a magnetic 
field is given by Lehnert in reference 7. It is found that a shallow 
layer of mercury over a copper disk with two concentric copper rings is 
destabilized by application of a vertical magnetic field. The rotation 
of the inner copper ring produces a shear layer in the mercury which is 
intensified by the magnetic field to the extent that an eddy- type flow 
results. It is pointed out by Lehnert that a generalization concerning 
the effect of a magnetic field on the flow field cannot then be made, 
and each situation must be studied to find out if the beginning of ampli- 
fication of a disturbance is actually delayed to a higher Reynolds number 
by the magnetic field. 

The flow of an incompressible electrically conducting fluid over a 
semi-infinite flat plate in the presence of a magnetic field perpendicular 
to the surface of the plate was studied in reference 8. The effect of 
the magnetic field on the stability of the flow has not as yet been 
studied for the case when the magnetic lines of force are perpendicular 
to or alined with the stream direction. It is the intent of this paper 
firstly to present an analysis of the stabilizing effect brought about 
by a coplanar magnetic field acting on an electrically conducting fluid 
flowing over a semi-infinite flat plate. The analysis is restricted to 
infinitesimal sinusoidal disturbances of the Tollmien-Schlichting type. 

In the course of the investigation it is necessary to evaluate a large 
portion of the numerical work for the corresponding two-dimensional chan- 
nel problem. Since the method of analysis is slightly different from 
that of reference 4, these results are presented. Secondly, the effect 
of a transverse magnetic field is considered. As was found for the chan- 
nel (Lock, ref. 5), the change in the critical Reynolds number,for the 
flat plate is controlled primarily by the change in the velocity profile 
brought about by the interaction of the fluid and magnetic field. The 
velocity profile shapes which are considered are taken from the two 
simplest cases analyzed in reference 8. The first case assumes that 
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the transverse magnetic field is fixed relative to the plate and the 
second that it is fixed relative to the fluid far from the plate. 

The method of analysis which is used is patterned after the procedure 
developed and described hy C. C. Lin in references 9, 10, and 11. A his- 
tory of the development and of the various physical problems which have 
been studied is given in a monograph by Lin in reference 12. A brief 
outline of the method is given in the introduction to the present analy- 
sis. The neutral stability curves are presented for several values of 
the magnetic parameter. 


SYMBOLS 


a 

B 

c 

F(z) 

m 


P 



imposed magnetic induction 
wave speed of disturbance 
Tietjen's function (see eq. (19)) 

crB 2 

magnetic parameter, , per unit length 


pressure 


lo> 0.1 A 

• • •/ 

R 


R x* 


u 

U 

U 


v 




in viscid perturbation amplitude functions (see eq. (8)) 


Reynolds number based on boundary-layer thickness, 


SU„ 


Reynolds number based on distance from leading edge of 

X*Uoo 

flat plate , — - — 
x component of velocity 

_U_ 

Uoo 

velocity in the stream direction of the flow field to be 
perturbed 

y component of velocity 
rectangular coordinates 
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a 

5 

€ 

♦ 

9 

v 

n 

P 

a 

j(o)> x (i) A 


00 


O 

d 


* 


y 0 (U 0 'aR) 1/3 

wave number of disturbance 

boundary-layer thickness, 6 
1 

(aR ) 1/3 

perturbation stream function 
amplitude function 
kinematic viscosity 

y-y Q 

e 

density of fluid 
electrical conductivity 

viscous perturbation stream functions (see eq. (l6)) 
Subscripts 

edge of boundary layer, or free stream 
critical layer where U = c 
disturbance 



Superscripts 


vector 

derivative with respect to y 
dimensional quantities 
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ANALYSIS 


The present state of stability theory requires that a number of 
simplifications be made in the analysis so that the method can be applied 
to physical situations without a prohibitive amount of labor. The method 
developed by C. C. Lin (refs. 9> 10, and 11) is a compromise between 
accuracy and effort required to analyze a given flow field. The present 
analysis is therefore patterned after it. 


Resume of Steps in Analysis 


The desired result is a stability diagram of exciting wave number a 
and Reynolds number R. At the beginning, the undisturbed steady- state 
solution to the magnetohydrodynamic flow problem being considered is 
assumed to be known. This information together with the equations of 
motion, the continuity equation. Maxwell 1 s equations. Ohm’s law for a 
moving fluid, the electromotive force relation, the wave nature of the 
disturbance, plus various approximations go to make up a complex fourth- 
order ordinary differential equation for the amplitude function cp. The 
various steps will now be explained. Sketch (a) was designed to orient 
the reader in the subsequent analysis which, in view of its well 
established nature, is discussed only briefly. 

The flow field is at some time assumed to be a steady two-dimensional 
stream of incompressible electrically conducting fluid. A two-dimensional 
infinitesimal sinusoidal disturbance of a given wave number a is then 
impressed on the fluid to test for the stability of the stream. A sinus- 
oidal disturbance is chosen because many disturbances which are likely 
to occur in nature can be Fourier analyzed and thereby reduced to a sum 
of sinusoidal disturbances. The magnitude of the disturbance is assumed 
to be vanishingly small or infinitesimal so that the analysis may be 
simplified by retaining only those terms which are linear in a disturbance 
or perturbation quantity. The wave nature of the disturbance is intro- 
duced by the disturbance stream function 


f = 5U 00 cpe ia *( x *- c ^ t - ic i n ) 


where a* is the wave number, c* is the velocity of the wave in the 
stream direction, and Cj_* is the rate of growth of the wave amplitude. 
The disturbance velocities are then given by u* = cty/dy* and 
v* = -(8i(r/8x*). The starred quantities have physical dimensions, whereas 
the unstarred counterparts have been made dimensionless by dividing by 
the free-stream velocity Uoo or by the boundary- layer thickness 6 as 
the case may be. It is assumed that the disturbance velocity and 
magnetic-field components are characterized by this exponential and 
depend on it to a first power. The object of the analysis is to find 


6 


NACA TN 4282 


the conditions -when the wave will just "begin to grow 1 (i.e., will "be 
neither damped nor amplified but neutral) , the factor cp* of the expo- 
nential is set equal to zero and the exponential reduces to ei a '*( x * -c *'t) . 
This function describes the propagation of the wave in the stream direc- 
tion for a given station x* as a function of time. The velocity of 
propagation for a given disturbance is independent of the distance along 
and perpendicular to the plate. It remains to find the circumstances 
■under which the wave amplitude neither grows nor diminishes but is neu- 
tral. It is found that the neutral point of wave growth occurs when the 
wave speed c* is equal to the local velocity U of the fluid. The 
region in the fluid where this happens is referred to as the "critical 
layer" and the distance from the wall as y Q *. 

In the actual flow problem one knows that the disturbance may be 
of either the two- or three-dimensional type. It has, however, been 
shown by Squire (ref. 14) that if the flow field is unstable to a three- 
dimensional disturbance it will be unstable to a two-dimensional disturb- 
ance at a lower Reynolds number. The extension of this proof to the type 
of magnetohydrodynamic problems being considered here is made by Michael 
in reference 15. Only two-dimensional disturbances will then be considered 
because they are the most unstable . 

When the information just described is introduced into the equations 
relevant to the problem, a complex ordinary fourth-order differential 
equation is obtained for the amplitude function cp (sketch (a)). 

It is complex because imaginary quantities are introduced by the expo- 
nential used to describe the perturbations. The terms which contain 
products or squares of the disturbance quantities are discarded. It is 
also assumed that the station in question is far enough downstream so 
that the variables are not changing in the free -stream direction. 

Even though a number of simplifying assumptions are made, the form 
of the differential equation is such that a simple solution has not yet 
been found. It is necessary then to find four linearly independent solu- 
tions by reducing the complete differential equation to two simpler dif- 
ferential equations by a power series expansion in l/R and e = (l/aR ) 1/3 
as indicated in sketch (a) . The zero-order terms in l/R are the only 
ones retained. The resulting differential equation is sometimes referred 
to as the inviscid form of the differential equation because all terms 
involving viscosity have been dropped. Proceeding down to the next step 


1 The stability curves corresponding to a number of growth rates, 
cp > 0 , have been computed by S. F. Shen (ref. 13 ) for flat plate and 
channel flow using an extension of Lin's method. 
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Complex fourth- order ordinary differential equation 
for 

<p(y,a,c,ir,m6,R) 


(l/ R ) 



1/R 




'I l/R*' / 3 e / 2 6 “d/aR) 

i Vi. v i 

(Discard terms involving viscosity) 


1/3 


Complex second-order 
differential equation ■ 
inviscid functions: 
cp(y,cc,c,U,m&) 




Complex fourth-order 
differential equation - 
viscous functions: 
X[(y-y 0 )(aR) 1/3 j Uo*] 



Cp 2 

X 1 = 

(y-y 0 )(a£) l/3 

X 2 = 1 

X 3 = 93 

*4 = 94 



Violates 
cp(S) = 0 
because 
tends to 
x) with y 


a,R combinations for neutral 
stability obtained from constants 
in 

cp = a^cpi + a 2 92 + &393 
with constants chosen so that 
cp(0) = cp 1 (0) = cp* (1) = 0 



Sketch (a) 


in sketch (a) , the first two linearly independent solutions cp]_ and cp 2 
are found "by introducing another series which consists of positive powers 
of the wave number a and whose coefficients depend on the wave speed c 
and velocity U in the flow field. Once again, only the first few terms 
in a are retained. 
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Hie boxes on the right of sketch (a) indicate that the solution to 
the first-order term in the set of differential equations -which results 
from the expansion in e = (l/aR) 1 ^ 3 is the only one which is found.. 

As pointed out by Lin, higher order functions could be found by quadra- 
tures but in most cases sufficient accuracy is obtained by considering 
only the first-order term. This differential equation has four linearly 
independent solutions which can be used. Two of these solutions, X x 
and X 2 , are discarded because they are too simple in form for curved 
velocity distributions. An examination of X 4 or cp 4 shows that it 
increases without limit with y and thereby violates the boundary con- 
dition that disturbances must die out as y approaches infinity. For 
this reason it is not used in the problems treated in this paper. The 
function cp 3 is generally used in the form known as Tietjen's function 
(sketch (a)). 

The three remaining linearly independent solutions are then combined 
in such a way that the boundary conditions are satisfied. The disturbance 
velocities will vanish at the wall and edge of the boundary layer when 

cp = a^cpi + a 2 cp 2 + a 393 

with the constants a!, a 2 , and a 3 chosen so that 

cp(0) = cp' (0) = cp' (1) = 0 

This is possible only for a certain combination of a and R when the 
magnetic parameter m5, velocity distribution U, and wave speed c have 
been specified. The end result from several such computations is a graph 
of the wave number a versus the Reynolds number R for various values 
of the magnetic parameter mS. Since these curves denote the values 
of a and R for neutral stability of the wave, a combination of a and R, 
which lies on the side of the curve denoted as unstable, warns that the 
amplitude of the disturbance will grow under those conditions. In the 
stable region the wave is damped. 

The number of approximations which are made might cause one to doubt 
the accuracy of the end results. Estimates made by Lin in reference 11 
indicate that the stability curves should not be in error by much more 
than a few percent and are therefore accurate enough for most engineering 
purposes. 

The analysis of the problems being considered in this paper is 
presented in the following sections. Since the method is well defined 
in references 9, 10, and 11, only the essential parts of the analysis 
are presented. 
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Coplanar Magnetic Field 


Differential equation . - The differential equation for the function cp 
will now he derived for the magnetic -field lines alined with the stream 
direction. The result is general enough that it can be applied to the 
flow in channels and over flat plates. Maxwell's equations for the 
incompressible -flow problems being considered are 


Div E = 0 
Div H = 0 


Curl H 
Curl E 







( 1 ) 


where E, H, j , and p are the electric field intensity, magnetic inten- 
sity, electric current density, and magnetic permeability, respectively. 
Ohm's law for a moving fluid is 

j = cr(E + U x B) (2) 

"where B = \iE and U is the local velocity vector. The equation of 
continuity is 


Div U = 0 (3) 

The Navier-Stokes equation modified to include the electromotive 
force term (so-called Lorentz force) arising from the relative motion 
between the fluid and magnetic field is 

+ (U • grad)U - - (j X B) + - grad p = vV a U (1*) 

dt P P 

where the excess charge density and applied electric field are assumed 
to be zero. 

The relation between the input wave number, a, of the disturbance 
and the Reynolds number, R, of the flow at which the amplitude of the 
disturbance neither increases nor decreases (neutral) will be found by 
introducing the quantities, 
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u* = U + 

v* = Vd *(y*) e ia.*(x»-c*t) 

B X = B + t(y*)e ia *( x5t - c * t ) 
By = 

F* - P* + 


\ 




(5) 


■where c* = c r * + xcj_*. When the disturbance is classified as neutral 
(neither amplified, nor damped.) c-j_* is zero. Since the problem will be 
to find only the neutral disturbance curves, the quantity c* will here- 
after be used to denote only the real part c r *, that is, the wave speed 
of the disturbance. The quantity a* is the wave number of the 
disturbance. 


It will be assumed that: 

1. The location of the instability is far enough downstream of the 
entrance to the channel or leading edge of the plate that the velocity 
normal to the boundary is negligible in comparison with the velocity 
parallel to the boundary. 

2. The fluid is of uniform density and conductivity, and the 
applied magnetic field, B, is uniform throughout the flow field. 

3. The boundaries are perfect conductors in order to complete the 
circuit for electric currents in the fluid. 

4. Terms which contain products or squares of the disturbance 
quantities are negligible. 

5. The disturbances are neutrally stable at values of the Reynolds 
number high enough so that a series in (l/aR) 1/3 converges rapidly. 

Following the method used by Stuart in reference 4, equations (l) 
through (5) may be combined and simplified using the foregoing assump- 
tions to yield a complex ordinary differential equation for the dimen- 
sionless amplitude function cp. 

(U-c) (q^'-a^cp) - cpU" + imSacp = -i- (cp ,m -2a 2 cp"+a 4 cp) (6) 

iaR 

where m = aB 2 /pU 00 , R = bU^/v, and U denotes the local velocity divided 
by the velocity at the edge of the boundary layer, Uo,. The symbols a 
and c in equation (6) denote the dimensionless form of the wave num- 
ber a* and wave speed c*, respectively. The amplitude function cp 
is a function of y = y*/6. The primes denote differentiation with 
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respect to the distance y normal to the nearest hounding surface. 
Hereafter, only the dimensionless uns tarred quantities will he used in 
the analysis unless it is noted otherwise. 

The boundary conditions are, 


cp = cp 1 =0 

at 

y = 0 



rat 

center 

of channel, y = 1 

0 

11 

'fe- 

ll 

0- 

-< or 
^-at 

y = 00 

for the flat plate 


Four linearly independent solutions to equation (6) will now he 
found hy the technique explained in references 9, 10, and 11. The first 
two solutions, cpj_ and cp 2 , will he derived from a series expansion in l/R 
and are designated as the inviscid solutions. The two remaining solutions 
93 and cp4 result from a series expansion in e = (l/aR) 1/3 and are 
called the viscous solutions. 

Inviscid solutions .- If the terms involving l/aR in equation (6) 
are assumed small, the remaining terms constitute the differential 
equation which cpi and cp 2 must satisfy. 

(U-c) (cp’’-a 2 cp) - cpU" + imhacp = 0 (7) 

A solution to equation (7) is found hy the method of Heisenberg 
(see, e.g. , ref. 12). It is assumed that the solution is of the form 

cp = % + aq x + a 2 q2 + a 3 q 3 + . . . (8) 

When equation (8) is inserted into equation (7) and the terms containing 
the same power of a are equated, the following set of linear ordinary 
differential equations is found. 

(9a) 
(9h) 

n = 2, 3, • • • (9n) 

The two linearly independent solutions of equation (9a) are 




U-c ^ 


= 0 


q - 


% “ 


U" 

U-c 

U" 

U-c 


qi = - 


imS 

U-c 


qo 


im5 

q n - qn-2 “ ^7 9n-i> 
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Obi = U - c 


(10a) 


% 2 



(101) 


The two inviscid solutions for the function cp may then be written 


as 


cpx(y) = (u-c) 


1 + 


r y r yi P 

a 2 / - ■ ■ 2 - / (U-c) 2 dy 2 dyi + 

Jo (u-c) Jo 

^1 0^2 


ny -1 oji p./^ .. p ys 

a 4 / -g- / (U-c) 2 / / (U-c) 2 dy 4 dy3dy2Clyi + . . . 

Jo (U-c) J 0 Jo (U-c) Jo 

(»-0) {taSa * 

ny p^l p^2 n p^3 

a 2 / — / (U-c) 2 / -5 / (U-c)dy 4 dy3dy 2 dy 1 + . . 

Jo (U-c )“ Jo Jo (U-c) 2 J 0 


* * j + ' 


( 11 ) 


and 


<te(y) = (u-c) 


im 5 a 


xx? * 

XXfC^Cx? *•••]*•••}♦••• as) 


where only linear terms in mb have been retained. 

The integrals in equations (ll) and (12) may be changed to a more 
convenient form by the transformation employed by Lin in reference 11. 
At the wall, y = 0, and at the edge of the boundary layer, y = 1, the 
inviscid functions and their derivatives then become 


9i(0) 

<Pj.(0) 

cp 2 (0) 


«Pg(0) 


-C 

U'(0) 

0 

_ 1 

c 


(13) 
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where 


cp 1 (l) = — im&<x(l-c)p 2 + . 


l-a 2 H 


cp' (1) = — 

' 1-c 


a 2 H 1 -a 4 M 3 + 


l-a 2 H 2 


- im5a(p 1 +a 2 p 3 + . . . ) 


92(1) = (1-c) 


K, 


_l-a 2 H 2 


a 2 N 3 + im&a(q 3 -KxP 2 - • • •) 


> (u) 


„•(!) = _1_ | -a*K 1 (H 1 -a*M3+ . . ,) + 1 . ^ . 

2 1-c \ l-a 2 H 2 


im5a 


Ki(Pi-Kt 2 p 3 + . . .) - P 2 - • • 


J 


Pi = / (U-c)dy 
Jo 


P 2 = r — 2 f (u-cjdysdyx 
Jo (U-c) J 0 


(15a) 

(15b) 


P3 


r 1 o r yi i r 72 

= / (U-c) 2 / -2 / (U-c)dy 3 dy 2 dy! 

Jo Jo (U-c) Jo 


(15c) 


13 = 


11 1 


lo (U-c )‘ 


pyl pi 

I (U-c) / • 

o J y 2 


(U-c)' 


dy 3 dy 2 dy x 


(15a) 


Hi = /"\u-c) ; 

Jq 


ay 


(I5e) 


H 2 = 


- 1 r y v ,2 


i 0 (u-c) J 0 


(U-c) dy 2 dy x 


(I5f) 


Ms 


<i5g) 
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N 3 



1 1 
(U-c) 2 



yi 


(U-c) 2 



p dy 3 dy 2 dyi 

(U-c) 2 


(I5h) 


Ki = 



1 1 
(U-c ) 2 


dy 


(I5i) 



Sketch (h) 


The path of integration, according to 
reference 10, lies along the real y axis with 
an indentation along a semicircular path under 
the singular point, y = y Q (i.e., where U = c) 
as shown in sketch (b) . 

Viscous solutions .- The two remaining 
independent solutions, the so-called viscous 
solutions cp 3 and cp 4 , are found by introducing 
the small parameter e = (l/aR) 1/3 and the 
function X as 


en = y - y Q 


<p(y) = xU) + eX, . (tj) +.e 2 X (q) + • • • 

(°) U) (2) 

and 

U” 

U - c = u£(er)) + n ( eT )) 2 + . . . 


(16) 


The subscript o indicates that the quantity is to be evaluated at the 
point where U = c. If the equations (l6) are introduced into equation (6) 
and the terms containing the same power of e are equated, the following 
set of ordinary differential equations results. 


IT 

• v nit 

+ lX 

= 0 

(°) 

(o) 


TT 

(i) 

4. iY"" 

x (i) 

= UqX 



Ui’X" , + im5X" 

2 ° (o) (o) 


(17a) 

(1Tb) 


etc. 

The solutions to equation (17a) are the only ones in this series 
which are found. As pointed out in the introduction to the analysis 
two of these four linearly independent solutions are discarded on the 
grounds that they are trivial. It is also found that the function X 4 
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or cp 4 increases with y indefinitely and thereby violates the boundary 
condition that the disturbance velocities must die out at the edg<* of the 
boundary layer. The form of the solution required is then 


<te(o) 

<Pg(0) 




( 18 ) 


where 


Z = y 0 (u’oR) 


1/3 


F(z) = 


/■v 


d£ / £ 1/2 H 


1/3 


f dU 3/2 




(19) 


-z J £ l/2 H 1/3 


3 (i ^ } 


3/2 




and Hi/ 3 ( ) is a Hankel function of the first kind and of order 1/3 . 
The function F(z) is sometimes referred to as the Tiet jen 1 s function. 
The tabulated values of references 9 a* 11 ! 16 are plotted in figure 1. 
The viscous solution is not modified by the presence of the magnetic 
field to the order of accuracy of the analysis. 


The inviscid solutions, equations (13) and (l4) , together with 
equation (18) make it possible to find the change in the neutral 
disturbance curve caused by a coplanar magnetic field. 


Channel flow (parabolic velocity profile) .- The flow of a viscous 
fluid between parallel planes (Poiseuille flow - sketch (c)) gives rise 
to a parabolic velocity profile if the station in 
question is not near the entrance to the channel. 

The effect of a coplanar magnetic field on the 


__r=£r-!= 


growth of a two-dimensional disturbance has already 
been studied by Stuart in reference 4. The differ- 
ence between the analysis carried out here and in 
reference 4 lies in the larger number of terms Sketch (c) 

retained here for the inviscid solution cp^. and in 

the form of equation (22) which is used to find the proper a-R combina- 
tion. The end results of the two analyses should, however, be about the 
same. Since the integrals (15) must be evaluated for a parabolic veloc- 
ity profile in order to make application to the flat-plate flow field, 
only a small amount of additional effort is required to find the neutral 
disturbance curves for the channel. 
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The velocity distribution is written as 

U = 2y - y 2 (20) 

The integrals (15) can he evaluated in closed form for arbitrary values 
of the wave speed c. 


Pi = - - 


If 1 , l-a 2 3a 2 -l , l+a\ iri / 0 2 o 3 

p 2 = — ( — + Zn + - — — Zn + (3a 2 -2a 3 -l) 

6 \a 2 a 2 2a 3 l-a/ 12a 3 


(21a) 

(21b) 


P3 = 2«i£ ( 3 a 2 -2a 3 -l) + Mli ln a + 

^ 45 18 V 5 4a a 300 ) 


1 ~ 17a//4 - (^l-3a + — ^ Zn(l+a) 


2a 3 +3a 2 -l 

60a 3 


(l+a) : 


2Ua a 


180 


1+2a3 ~ 3a2 (l-a) 3 
60a 3 


211af + l+l 7 a/4 _ / + 8af) , n(l . a) 

L 180 5 V 3 / 


(21c) 




J ^ 1+a ^ - (2a-l) Zn(l+a) + (2a+l) Zn(l-a)+3(a 2 +2) 


2a 


2a 


L » Os) - - \ - ^ • Ii 2 v ilj! * <- 


i it 


16a 6 L 3 


3S1 zn(l-a 2 ) - y in 2a + 2a(2-5a g )_ + 2 ^ a 2_ 1^ ln (lt0 + ( 3 a 2 -2a 3 -l) 


4a 3 


(2ld) 


where a 2 = 1-c, and L 2 ( ) is the dilogarithmic integral. Numerical 
values for the relations (21) for several values of the parameter c are 
presented in table I. The functions L 2 are tabulated in references 17 
and 18. The remaining integrals in the group (15) are written and 
tabulated in reference 11. 

It remains now to combine the inviscid and viscous solutions so that 
the boundary conditions at the wall and at the edge of the boundary layer 
are satisfied. It is found that the wave number of an antisymmetric 
disturbance and the Reynolds number of the flow field must be chosen so 
that 
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F(z) 


<P 3 (°) 

y^(o) 


c<Pg(l) 


u 1 (o)cp^(l) + i cp^(l) 


( 22 ) 


With all of the individual functions known, an iteration scheme is 
employed to find the correct wave number and Reynolds number combination. 
A graphical method was used to find the intersection of the curves of 
the functions on the left and right sides of equation (22) for several 
values of the parameter c, whereas a numerical iteration scheme was 
used in reference 11. 

The neutral disturbance curves for several values of the magnetic 
parameter mS are shown in figure 2(a). Since the parameter m5a was 
held constant in the analysis of reference 4, a direct comparison with 
the neutral stability curves of that paper cannot be made. The critical 
Reynolds numbers found by the two analyses will be compared in the 
discussion. 


Flat-plate velocity profile.- When an incompressible viscous fluid 
flows past a semi- infinite flat plate of zero thickness, the velocity 
profile can be predicted theoretically and 
is generally referred to as the Blasius pro- 
file (sketch (d)). The neutral stability 

curve in the nonmagnetic case has been com- U« 
puted in references 9* 10, and H« The 
effect of the magnetic field on these 

results will now be found. Sketch (d) 




U — 

B 

— 


— 

_ — — " 






The integrals (15a) through (l5d), evaluated by the approximate 
method suggested in reference 11, are tabulated in table II for specific 
values of c. The real and imaginary parts of K x are computed by the 
relations given in reference 11 as 


K 


ir = 


cU» (0) 


+ 0.1465 + 1.2467c + 1.045c 2 + 2.039c 3 + 4.078c 4 + 


2.423c 5 + . . . + .2 ( c 2 + -2± c 5 + . . . 


8 
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8 


in 2l§l£ + ijt 


- 


U" 

-rt 




(23) 

(24) 


The expression which determines the proper values of wave number a 
and the Reynolds number R for the flat-plate problem is 
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F(z) 


c[cpg(l) + cup g (l)] 

y Q |u' (0) [cp^(l) + cwp 2 (l)] + | [q>i(!) + 0^(1)]} 


(25) 


The neutral stability curves for several values of the magnetic 
parameter mb are shown in figure 2(b) . 


Transverse Magnetic Field 


The change in the boundary- layer velocity profile for flow over a 
flat plate in the presence of a transverse magnetic field was found in 
reference 8. It was found that the skin friction and heat transfer are 
reduced if the magnetic field is fixed relative to the plate (sketch (e)) 

and increased if it is fixed relative to the 
fluid outside of the boundary layer 
(sketch (f)). The possibility exists, however, 
that the magnetically induced velocity profile 
may be more or less stable to transition to 
turbulent flow. An estimate of the change in 
the stability of an infinitesimal sinusoidal 
disturbance induced by the transverse magnetic 
field will now be found. 


The differential equation for the 
disturbance stream function is found by the 
technique used by Lock in reference 5 which 
is to combine equations (l) through (5) and 
then simplify the result by applying the five 
assumptions outlined in the analysis of the coplanar field. The differ- 
ential equation for the perturbation stream function is then 

(U-c) (cp"-a 2 cp) - U"cp = cp" + -L ( cp"" -2a 2 cp"+a 4 cp) (26) 

a* iaR 



|M 11 

U« 

-M 1 

ill 

1 |U 

-l—u 

Ml 

Ij-J-H 

11 

H'l't 

M 


Sketch (f) 


It is shown by Lock in reference 5 that the forms of the inviscid 
and viscous solutions are not affected to the order of the analysis by 
the additional magnetic term in equation (26) . In other words, the 
change in the velocity profile caused by the transverse magnetic field 
dominates the stabilizing action of the magnetic field. The neutral 
stability curves for several values of the magnetic parameter mx are 
found by the method outlined in the appendix of reference 11 . The invis- 
cid solutions are found by using the numerical data in tables I and II 
of reference 8 to determine the velocity profiles at mx = 0.05 and 0.10. 
The numerical results for the integrals (l 5 e) through (I 5 h) are tabulated 
in tables III and IV. 
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The real part of the integral (l5i) is evaluated by expanding in a 
series about the critical point y = y 0 where U = c. The result which 
was used in the computations for the transverse magnetic field is 


K-, 


0.4 

y 0 (o.i.-y 0 )(uy 2 


IT" 

u o 


(u’) : 


in 


0 - 4 -y o 


+ 


0.3 


(Up ) 2 

( u ') 4 


( u ") 3 

— (0.04-0. 2y„) + 

( U ^) 5 


5 ( u p) 4 ( 0.064 



- 0 . l6y 0 + 0.4y o 


4(l-c) V 0.75-c \/i-c 


Zn 



0.1 

(l-c) ; 


(27) 


The imaginary part K 1± is evaluated by use of equation (24) . The 
velocity U in the integrals (I5e) through (I5h) , (24) , and (27) is 
referred to the velocity at the edge of the boundary layer a L the partic- 
ular station being considered. When the magnetic field is fixed relative 
to the plate the undisturbed stream velocity and the velocity at the edge 
of the boundary layer are not the same. 

The neutral disturbance curves are shown in figures 3(a) and 3(b). 


DISCUSSION 


The neutral stability curves shown in figures 2 and 3 indicate that 
the presence of a magnetic field may stabilize or destabilize the flow 
of an incompressible, electrically conducting fluid. It is seen from 
these results that the flow over a flat plate is stabilized by either a 
coplanar magnetic field or by a transverse magnetic field fixed relative 
to the fluid, but a transverse magnetic field fixed relative to the plate 
is generally destabilizing. The portion near the top of the mx = 0.1 
curve in figure 3(a) indicates an opposite trend for a small range in 
wave number. As pointed out in the introduction, another example of 
flow instability caused by a magnetic field is presented by Lehnert in 
reference 7» 

A given flow field will probably contain disturbances covering a 
wide range of wave number due to imperfections in the walls and entrance 
to the flow field. A conservative value for the critical Reynolds number 
is then the lowest value at which it is first possible for any of the 
waves to be amplified. The critical Reynolds numbers for the flow prob- 
lems considered in references 4 and 5 and for the coplanar magnetic- field 
cases studied in this paper are shown in figure 4 as a function of the 
magnetic parameter m8. The results for the transverse magnetic field 
as a function of mx are also shown in figure 4. It is seen that the 
results of Stuart in reference 4 are in essential agreement with the 
present analysis. The difference between the results is attributable 
to the smaller number of terms retained in the analysis of reference 4 
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for the inviscid solution. The results for a laminar mixing region 
obtained by Curie in reference 6 are not shown in figure 4 because the 
Reynolds numbers are too small for the scale of the graph. 

It is quite evident from figure 4 that a magnetic field is more 
effective when applied to channel flow than to flat-plate flow. In 
particular, the transverse magnetic field is so effective in stabilizing 
the flow in a channel that the curve is a vertical line to the scale of 
the graph. 

When the magnetic field is coplanar, the large difference in the 
shape of the critical Reynolds number curves for the channel and flat- 
plate flow fields is attributable to the infinite extent of the flow 
field above the flat plate. As is shown by Lin in reference 9, the 
asymptotic form of the disturbance stream function as the distance y* 
approaches and exceeds the boundary- layer thickness, 8, introduces addi- 
tional terms in the equation determining the neutral stability curves. 
This is obvious when equations (22) and (25) are compared. These addi- 
tional terms de-emphasize the terms involving the magnetic parameter and 
result in a much smaller stabilizing effect for the flat plate than for 
the channel flow. 

The magnetic parameter and the Reynolds number for the flow over a 
flat plate at which an infinitesimal disturbance will grow (figs. 2(b) 
and 3) are based on the boundary- layer thickness 5 taken 2 as 
6/ n/u to /vx*, where U/Uoo ~ 0.999* The distance along the plate from the 
leading edge is then related to the boundary- layer thickness by the 
relationship 


6x* 


where, R x * 


UoqX* / v . Therefore , 


and 


m8 


6mx* 

J Rx* 


R = 6 s/Rx* 


2 Stand ard tex ts on boundary -layer theory usually define the thickness 
as 8 = 5/ "VUco/vx*, where U/U^ ~ 0. 99* As explained in reference 11, 
more accuracy is achieved by defining a thicker boundary layer to a 
evaluate the inviscid integrals. 
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It was found in reference 8 that a magnetic field perpendicular to 
a flat plate changes the velocity profile in the boundary layer. Even 
a small magnetic field fixed relative to the plate will cause an inflec- 
tion point 3 in the velocity profile near the surface. As is shown in 
figures 3(a) and 4, this causes the flow to he less stable with a mag- 
netic field than in the nonmagnetic field case. The results in fig- 
ures 3(b) and 4 indicate that a magnetic field fixed relative to the 
fluid far from the plate changes the velocity profile to a shape which 
is more stable. The results of reference 8 indicate that the skin fric- 
tion and heat transfer are reduced in the former and increased in the 
latter case. Care must then be exercised if one attempts to reduce either 
the skin friction or heat transfer by imposing a magnetic field across 
(perpendicular to) the flow field and not in relative motion with the 
plate, because the laminar flow is destabilized by this technique. Like- 
wise, the increase in the skin friction and heat transfer brought about 
by a transverse magnetic field sweeping past the plate at the velocity 
of the free stream would eventually experience a moderate compensating 
effect in the form of increased stability of the laminar stream. 

The results of this paper, in conjunction with that of reference 8, 
point out the fact that it is not certain whether the skin friction and 
heat transfer are lowered or raised by using a transverse magnetic field 
to alter the flow over a flat plate. The magnetic field alters the 
velocity profile and changes the rate of growth of small disturbances so 
that the two effects tend to compensate each other. Individual situations 
must then be considered separately to determine whether an advantage can 
be achieved. 


CONCLUSIONS 


The analysis carried out in this report for the flow over a flat 
plate indicates the effect of a magnetic field on the stability of a 
disturbance of the Tollmien-Schlichting type. In particular it is found 
that : 


1. The flow is stabilized by a coplanar magnetic field. The 
increase in the critical Reynolds number is small compared with the 

increase achieved in a channel with a coplanar or transverse magnetic 
field. 


2. A transverse magnetic field fixed relative to the flat plate 
changes the velocity profile to an inherently unstable shape which lowers 
the critical Reynolds number. 


It is noted in figure 3(a) that the maximum value of the wave number 
first increases and then decreases with increasing mx. This is caused 
by the rapid change in the curvature of the velocity profile with mx. 
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3 . A transverse magnetic field fixed relative to the fluid far from 
the plate changes the velocity profile in the boundary layer to a shape 
which is more stable and thereby raises the critical Reynolds number. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 1, 1958 
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TABLE I.- COEFFICIENTS FOR INVISCID SOLUTIONS; PARABOLIC 
VELOCITY PROFILE (COPLANAR MAGNETIC FIELD) 


c 

Pi 

P2 

P3 

03 

0 

0.6667 

0.3977 

0.09281 

00 

.05 

.6167 

.4101 

.0858 

0.6765 

.1 

.5667 

.4224 

.07900 

.6529 

.15 

.5167 

.4351 

.07247 

.6675 

.2 

.4667 

.4480 

.06617 

.7035 

.25 

.4167 

.4615 

.06010 

• 7587 

.3 

.3667 

.4756 

.05424 

.8355 

.35 

.3167 

.4905 

.04855 

.9401 

.4 

.2667 

.5062 

.04302 

1.0828 

.45 

.2167 

.5230 

.03762 

1.2818 


TABLE II.- COEFFICIENTS FOR INVISCID SOLUTIONS; BLASIUS 
PROFILE (COPLANAR MAGNETIC FIELD) 


C 

Pi 

P2 

P3 

03 

0 

0.7133 

0.3817 

0.1272 

00 

.05 

.6633 

• 3941 

.1202 

.5603 

.1 

.6133 

.4o64 

.1134 

.5367 

.15 

.5633 

.4191 

.1068 

.5512 

.2 

.5133 

.4321 

.10053 

.5872 

.25 

.4633 

.4455 

.09546 

.6424 

.3 

.4133 

.4597 

.08860 

.7192 

.35 

.3633 

.4745 

. 08261 

.8238 

.4 

.3133 

.4902 

.07739 

.9665 

.45 

.2633 

.5070 

.07198 

1.1656 
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TABLE III.- COEFFICIENTS FOR INVISCID SOLUTIONS; TRANSVERSE MAGNETIC 

FIELD FIXED RELATIVE TO PLATE 


mx 

c 

Hi 

h 2 

M s 

Ns 

Ki 

x r 

•H 

H 

Uo 


0.05 

0 

0.6022 

0.2337 

0.07907 

0.1700 

-OO 

0 

1.8273 

0 


.05 

• 5372 

.2225 

.06916 

.1806 

-10.87 

-.0310 

1.8284 

.0274 


.1 

.4772 

.2107 

.06052 

.1929 

-5.327 

-.0338 

1.8302 

.0547 


.15 

.4221 

.1983 

.05310 

.2044 

-3-414 

-.0102 

1.8312 

.0820 


.2 

• 3721 

.1854 

.04682 

.2158 

-2.394 

.0397 

1.8303 

.1093 


.25 

.3271 

.1718 

.04162 

.2266 

-1.721 

.1159 

1.8263 

.1366 


.3 

.2870 

.1576 

.03741 

.2385 

-1.209 

.2199 

1.8172 

.1641 


.35 

.2520 

.1428 

.03409 

.2466 

-.7729 

.3553 

1.8009 

.1917 


.4 

.2220 

.1274 

.03151 

.2472 

-.2995 

.5209 

1.7844 

.2197 


.45 

.1969 

.1112 

.02956 

.2389 

.2253 

.7250 

1.7509 

.2477 

.10 

0 

• 5733 

.2223 

.06975 

.1732 

-00 

0 

1.6473 

0 


.05 

.5089 

.2111 

.05984 

.1838 

-12.28 

-.0914 

1.6498 

.0303 


.1 

.4495 

• 1993 

.05120 

.1962 

-6.156 

-.1394 

1.6556 

.0606 


.15 

.3951 

.1869 

.04378 

.2076 

-4.051 

-.1456 

1.6625 

.0907 


.2 

.3456 

.1740 

.03750 

.2190 

-2.939 

-.1132 

1.6686 

.1207 


.25 

.3012 

.1604 

.03230 

.2298 

-2.228 

-.0447 

1.6724 

.1507 


.3 

.2618 

.1462 

. 02809 

.2417 

-1.717 

.0580 

1.6719 

.1806 


.35 

.2274 

.1314 

.02477 

.2498 

-1.323 

.1950 

1.6657 

.2105 


.4 

.1980 

.1160 

.02219 

.2504 

-1.0126 

.3726 

1.6520 

.2408 


.45 

.1736 

.0998 

.02024 

.2421 

-•7831 

• 5959 

1.6310 

.2718 
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TABLE IV.- COEFFICIENTS FOR INVISCID SOLUTIONS; TRANSVERSE MAGNETIC 

FIELD FIXED RELATIVE TO FLUID 


mx 

C 

Hi 

H 2 

% 

w 3 

Ki 

x r 

Kxi 

% 


0.05 

0 

0.6458 

0.2536 

0.09507 

0.1704 

-00 

0.4537 

2.3245 

0 


.05 

.5753 

.2424 

.08516 

.1810 

-7.937 

.4621 

2.2862 

.0218 


.1 

• 5099 

.2306 

.07652 

.1934 

-3.699 

.4748 

2.2481 

.0438 


.15 

.4494 

.2182 

.06910 

.2048 

-2.248 

.4936 

2.2101 

.0662 


.2 

• 3940 

.2052 

.06282 

.2162 

-1.4584 

.5319 

2.1707 

.0891 


.25 

• 3435 

.1917 

.05762 

.2270 

-.9174 

.5874 

2.1293 

.1124 


• 3 

.2981 

.1775 

.05341 

.2389 

-.4697 

.6642 

2.0856 

.1360 


• 35 

.2576 

.1627 

.05009 

.2470 

-.0392 

.7692 

2.0370 

.1603 


.4 

.2222 

.1472 

.04751 

.2476 

.4347 

.9107 

1.9826 

.1852 


• 45 

.1917 

.1310 

.04556 

.2393 

1.0287 

1.0983 

1.9219 

.2108 

.10 

0 

.6585 

.2604 

.10067 

.1725 

-00 

.6430 

2.6073 

0 


.05 

.5768 

.2492 

.09076 

.1831 

-6.648 

. 6606 

2.5390 

.0196 


.1 

.5202 

.2374 

.08212 

.1954 

-2.966 

.6869 

2.4712 

.0395 


.15 

.4586 

.2250 

.07470 

.2068 

-1.6996 

.7194 

2.4o45 

.0600 


.2 

.4019 

.2121 

.06842 

.2183 

-1.0201 

.7611 

2.3380 

.0812 


.25 

.3503 

.1985 

.06322 

.2291 

-.5522 

.8189 

2.2712 

.1028 


.3 

.3036 

.1843 

.05901 

.2410 

-.1525 

.8923 

2.2029 

.1252 


.35 

.2620 

.1695 

.05569 

.2491 

.2445 

.9911 

2.1320 

.1483 


.4 

.2253 

.1541 

.05311 

.2497 

.7030 

1.1244 

2.0584 

.1721 


.45 

.1937 

.1379 

.05116 

.2414 

1.3022 

1.3015 

1.9803 

.1968 
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Figure 1.- Tietjen’s function, F(z) = F r (z) + iFj_(z) where z = y 0 (U 0 T aR) ' . 
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(a) Channel flow; parabolic velocity profile. 

Figure 2 .- Regions wherein an infinitesimal sinusoidal disturbance is amplified or damped in the 

presence of a coplanar magnetic field. 
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(b) Flat plate; Blasius boundary- layer profile. 


Figure 2.- Concluded. 
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(a) Magnetic field fixed relative to plate. 

Figure 3.- Regions wherein an infinitesimal disturbance is amplified or damped for flow over a 
semi-infinite flat plate in the presence of a transverse magnetic field. 
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(b) Magnetic field fixed relative to fluid far from the plate. 
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Figure 3*- Concluded. 
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Figure 4.- Critical Reynolds number as a function of the magnetic parameter. 
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